Highly-conserved regulatory activity of the ANR family in the virulence of diarrheagenic bacteria through interaction with master and global regulators

ANR (AraC negative regulators) are a novel class of small regulatory proteins commonly found in enteric pathogens. Aar (AggR-activated regulator), the best-characterized member of the ANR family, regulates the master transcriptional regulator of virulence AggR and the global regulator HNS in enteroaggregative Escherichia coli (EAEC) by protein–protein interactions. On the other hand, Rnr (RegA-negative regulator) is an ANR homolog identified in attaching and effacing (AE) pathogens, including Citrobacter rodentium and enteropathogenic Escherichia coli (EPEC), sharing only 25% identity with Aar. We previously found that C. rodentium lacking Rnr exhibits prolonged shedding and increased gut colonization in mice compared to the parental strain. To gain mechanistic insights into this phenomenon, we characterized the regulatory role of Rnr in the virulence of prototype EPEC strain E2348/69 by genetic, biochemical, and human organoid-based approaches. Accordingly, RNA-seq analysis revealed more than 500 genes differentially regulated by Rnr, including the type-3 secretion system (T3SS). The abundance of EspA and EspB in whole cells and bacterial supernatants confirmed the negative regulatory activity of Rnr on T3SS effectors. We found that besides HNS and Ler, twenty-six other transcriptional regulators were also under Rnr control. Most importantly, the deletion of aar in EAEC or rnr in EPEC increases the adherence of these pathogens to human intestinal organoids. In contrast, the overexpression of ANR drastically reduces bacterial adherence and the formation of AE lesions in the intestine. Our study suggests a conserved regulatory mechanism and a central role of ANR in modulating intestinal colonization by these enteropathogens despite the fact that EAEC and EPEC evolved with utterly different virulence programs.


Characterization of the Rnr regulon in EPEC.
Since enteric pathogens use distinct mechanisms to colonize and invade their hosts, and the fact that the amino acidic identity between ANR members differs significantly among pathogens (Fig. 1B), it is uncertain whether ANR accomplishes the same regulatory function in different enteric pathogens. We sought to determine this gap in knowledge by dissecting the regulon and biological role of Rnr in attaching and effacing pathogens, which shares only ~ 25% of amino acid identity with the archetype Aar. Accordingly, the rnr gene was deleted in the prototype EPEC O127:H6 strain E2348/69, and RNA-seq determined its transcriptome. For these experiments, the wild-type (WT) EPEC, its isogenic EPECrnr mutant, and EPECrnr complemented with rnr in-trans [EPECrnr(pRnr)] were grown in DMEM-high glucose for six hours to activate the expression of Rnr. Subsequently, total RNA from all strains was extracted and processed for cDNA synthesis, library construction, and DNA sequencing by CD genomics (NY, USA), as indicated in the material and methods. Bioinformatics analysis revealed approximately 500 genes that were differentially expressed (DEGs) in the Rnr regulon (± 1.5 fold, P < 0.05) (Fig. 2C-F). The majority of Rnr-regulated genes were located in the chromosome of EPEC E2348/69 ( Fig. 2A) and associated with six major functional categories: genes involved in metabolism (46%), protein transport (11%), regulation (5%), ribosomal activities (7%), virulence (including bacterial adherence and motility) (3%), and other functions (28%) (Fig. 2B,G-L and Supplementary Fig. S1). 8% of Rnr-regulated genes were encoded in the LEE pathogenicity island located in the chromosome and 6% in the pMAR2 plasmid (GenBank FM180569.1) ( Fig. 2A).
Ler and HNS were identified in the Rnr-regulon with 26 other regulatory proteins ( Fig. 2G and Supplementary  Fig. S1). Most of these regulators belong to the AraC/XylS family, including PerA, EutR, MelR, AdiY, YdiP, and GadX. Interestingly, we previously found that Aar also regulates GadX in response to the acid environment in EAEC 30,31 . Other important regulators under Rnr control are BssG, FimG, and CsgD, associated with biofilms or bacterial adherence (Fig. 2G).
In addition, we observed that a large number of genes regulated by Rnr are involved in metabolism (~ 200 genes), including the Rut operon ( Supplementary Fig. S1), intricate in the degradation of exogenous pyrimidines as the sole nitrogen source, and the arginine succinyltransferase pathway which uses arginine as a source of carbon and nitrogen. Numerous genes of these operons are also under the NtrC control 32,33 . Rnr negatively regulates the locus of enterocyte effacement pathogenicity island (LEE-PAI). The LEE-PAI (~ 36 kb) is composed of 42 genes and seven operons encoding the T3SS (Fig. 3A).
Our transcriptomic data shows increased expression of LEE genes in EPECrnr, whose complementation in trans www.nature.com/scientificreports/ with pRnr plasmid restored the expression of genes to comparable wild-type levels (Fig. 3B). Among those, twenty genes encode for proteins that form the core of the T3SS were at least ~ 2-fold upregulated, including proteins of the basal body (EscC, EscD, EscJ); inner membrane machinery (EscV, EscR, EscS, EscT, EscU); needle tip and translocon (EspD, EspB, EspA) and the EscN ATPse. The intimin gene (eaeA) and its receptor tir were also upregulated (~ 5-fold), which are associated with EPEC intimate adherence to host cells (Fig. 3B). Several LEE regulators (Ler, GrlA and GrlR) located in the LEE-PAI were upregulated (~ 2 to 18 fold) in EPECrnr as compared to WT and complemented EPECrnr strain (Fig. 3B).
To validate our transcriptome dataset, nine genes were selected based on their relevance in the virulence of EPEC for qRT-PCR analysis. We analyzed genes encoding structural proteins of the T3SS apparatus (espA, espB, and espD), regulatory proteins (hns, ler, and perA), and proteins involved in adherence (tir, eaeA, and bfpA). For qRT-PCR experiments, WT EPEC, EPECrnr, and EPECrnr(pRnr) strains were grown in DMEM-high glucose for 2, 3, and 4 h, and total RNA was isolated and prepared for gene expression analysis. Our data showed that 8 out of 9 analyzed genes exhibited higher levels of expression in the EPECrnr mutant (~ 3-fold) compared to WT after 4 h of growth (middle log phase) (Fig. 4).
We previously found that Aar increases the expression of AAF genes in early EAEC growth stages by decreasing the expression of HNS, which acts as a repressor of AAF expression; however, when Aar is increased, it acts as a negative regulator of AAF by inactivating AggR, the positive AraC/XyS regulator of AAF 34 . Similar findings were observed in EPEC with bfpA and eaeA (Fig. 4G,I), which are regulated by PerA and HNS, respectively, and these, in turn, are regulated by Rnr (Fig. 4A,H).  www.nature.com/scientificreports/ We sought to determine if changes observed at the transcriptional level also correlated with changes at the protein level. Accordingly, EPEC derivatives were grown in DMEM until the late exponential growth phase (OD 600 ηm of 1.0), and the abundance of T3SS structural proteins (EspA and EspB) (

Rnr protein directly interacts with HNS and Ler proteins.
We previously showed that Aar interacts with HNS global repressor affecting its regulatory activity 30 . Therefore, we sought to determine if Rnr can interact with HNS and Ler, both members of the HNS family. The Bacterial-two hybrid (BACTH) system is broadly used to scrutinize protein interactions between regulatory proteins 35 , and we have successfully used this approach to examine interactions between Aar-AggR and Aar-HNS 30,31 . Thus, we used the BACTH system to investigate interactions between Rnr, HNS, and Ler. Accordingly, rnr, hns, and ler genes were fused to T25 and T18 fragments of the catalytic domain of Bordetella pertussis adenylate cyclase, expressed in plasmids pKNT25 and pUT18, respectively (Fig. 6A) 35 . The resulting plasmids were co-transformed in different combinations of pUT18 and pKNT25 derivatives into the reporter strain E. coli BTH101. Remarkably, we observed protein-protein interactions of Rnr with members of HNS family; HNS, and Ler in the BACTH system manifested by the appearance of a moderate to intense blue color on agar plates (Fig. 6C) and quantification of the β-galactosidase activity (Fig. 6B,D). Taken together, our findings suggest that Rnr is regulating gene expression of the LEE-PAI by direct interaction with Ler and HNS global regulators.

Rnr and Aar are interchangeable for regulating the T3SS in EPEC and fimbriae in EAEC.
Since the heterologous expression of Aar in EPEC is capable of downregulating the expression of T3SS in EPEC in our previous experiment (Fig. 5B,C), we sought to determine whether the heterologous expression of Rnr in EAEC042 downregulates the expression of the AggR-regulated AAF, the main virulence factor of EAEC associated with host-interactions. Accordingly, we analyzed whole-cell proteins from EAEC derivatives expressing Aar www.nature.com/scientificreports/ and Rnr by SDS-PAGE and Western blot (Fig. 7). As expected, we found that Rnr was able to downregulate the expression of the major AAF fimbria subunit, AafA, in EAEC (Fig. 7), suggesting that despite the low homology between Aar and Rnr, they may possess structural features that allow function conservation between distantly related ANR members.
Aar and Rnr negatively impact intestinal colonization in their respective pathogens. Despite the extensive molecular characterization of ANR in EAEC, its role in bacterial pathogenesis is not entirely understood, partly due to the lack of adequate animal models for E. coli pathogens. Since Aar downregulates the expression of AggR-regulated AAF in EAEC (Fig. 7) 29 , and Rnr downregulates PerA-regulated BfpA and genes associated with intimate adherence mediated by T3SS in EPEC (Fig. 4), we sought to determine the impact of ANR regulation in bacterial adherence and intestinal colonization. Human intestinal organoids have become the gold standard for studying host-pathogen interactions and have been successfully used to investigate essential features of EAEC and EPEC pathogenesis 36,37 . We, therefore, www.nature.com/scientificreports/ used this relevant intestinal model to examine the role of ANR in bacterial colonization. For these experiments, human intestinal colonoid monolayers were infected with parental 042, 042aar, and 042aar (pAar) at 37 °C for 6 h, and bacterial adherence was analyzed by confocal microscopy (Fig. 8A-P). We observed that the deletion of aar significantly increases bacterial colonization in human colonoids compared to the parental strain (Fig. 8C,G,K,O). In agreement with the negative role of Aar, complementation of 042aar with the pAar plasmid drastically reduced biofilm formation (Fig. 8D,H,L,P) as judged by the enumeration of bacterial cells on colonoids (Fig. 8R). Moreover, microscopic examination of bacterial biofilms revealed increased bacterial aggregation in colonoids infected with 042aar than parental 042 strain (Fig. 8N,O). We next determined whether Rnr impacts EPEC intestinal colonization and the formation of AE lesions in human intestinal organoids. Accordingly, human cell monolayers were infected with 10 6 CFU of EPEC E2348/69, EPECrnr, and EPECrnr(pRnr) for 6 h at 37 °C. Subsequently, infected cells and uninfected controls were analyzed for EPEC adherence and formation of AE lesions by confocal microscopy (Fig. 9). The confocal images were pixel-quantified as previously reported (Fig. 9M) 36 . We observed a more significant number of adhering EPECrnr strain on intestinal cell monolayers and which correlated with a greater number of AE lesions compared to the parental EPEC strain (Fig. 9A,B,D,E,G,H,M). Complementation of EPECrnr with the pRnr plasmid drastically reduced bacterial adherence (Fig. 9C,F,I) and the number of AE lesions on intestinal cell monolayers as judged by actin polymerization beneath the adherent bacteria (Fig. 9L). Although AE lesions were observed in cells infected with all EPEC strains, cells infected with EPECrnr(pRnr), which overexpresses Rnr, exhibit smaller actin pedestals than the WT or EPECrnr strains as judged by qualitative analysis of confocal images (Fig. 9J,K,L). Taken together, our data suggest a central role of ANR in modulating intestinal colonization by diarrheagenic pathogens (Figs. 8, 9).

Discussion
The ability of a pathogen to colonize a host and cause disease requires coordinated expression of genes that mediate nutrient acquisition, as well as genes involved in virulence. The mechanisms used by bacterial pathogens to regulate their fitness and virulence have been the subject of intense investigation for several decades. One such regulatory mechanism is exerted by the recently identified ANR family of small regulatory proteins whose mode of action is to interfere with the function of global regulators by protein-protein interactions, thereby modulating gene expression globally 27,[29][30][31] . In EAEC, the prototype ANR termed Aar directly interacts with AggR and HNS proteins, resulting in an impaired ability of these regulators to bind DNA and, consequently, affecting the expression of hundreds of genes associated with metabolism, stress-response, fitness, and virulence 27,30,31 .
In this study, we sought to determine whether the regulatory properties of ANR are conserved in AE pathogens, such as EPEC, which has a multifaceted virulence scheme and whose ANR (Rnr) shares only 25% identity with Aar of EAEC. Accordingly, we defined the Rnr regulon of classical EPEC strain E2348/69 by comparing the transcriptomes of WT EPEC, isogenic EPECrnr, and the complemented strain. We found that, like Aar in EAEC, Rnr regulates genome-scale gene expression in EPEC. More than 500 genes, including genes associated with fitness and virulence, were controlled by Rnr (Fig. 2). As expected, Rnr regulates the expression of HNS and multiple AraC/XylS regulators present in EPEC, including MelR, AdiY, YdiP, EutR, GadX, and PerA ( Fig. 2G and Supplementary Fig. S1), which are associated with carbon metabolism, stress-response, and virulence functions, respectively (Fig. 2B). Even more remarkable is that Rnr negatively regulates the expression of the entire LEE pathogenicity island, including genes required for expressing the T3SS and its positive LEE-encoded regulator (Ler) (Fig. 3, Supplementary Fig. S1). Ler belongs to the HNS family, but unlike other HNS members, Ler activates the transcription of LEE genes by counteracting HNS-mediated repression 38 . We found that the expression of genes encoded in LEE2, LEE3, and LEE5 under the HNS-Ler control was downregulated by Rnr as judged by RNAseq and qRT-PCR transcriptional analysis (Figs. 3 and 4). Moreover, we observed that the abundance of T3SS-structural proteins EspA and EspB (encoded in LEE4) was moderately lower in whole-cell and supernatant preps of the WT strain when compared to EPECrnr but drastically reduced in the complemented strain overexpressing Rnr (Fig. 5), highlighting the negative regulatory role of Rnr.
On the other hand, we confirmed the direct protein-protein interaction of Rnr with Ler and HNS in vivo using the BACTH system, which has been successfully used to test interactions between regulatory proteins, including Aar 27,30,31 . Since Rnr binds to both a silencer (HNS) and anti-silencer (Ler), and the fact that binding of Aar to HNS and AggR in EAEC hampers the function of these regulators, we have envisioned two possible mechanisms that may explain the negative regulatory role of Rnr on the expression of LEE genes: (1) direct Rnr interaction with Ler may inhibit its positive feedback activity on the LEE PAI and (2) since HNS function as a repressor of its transcription, binding of Rnr to HNS may inhibit the HNS silencing effect on its promoter, favoring HNS expression, and consequently, the repression of LEE genes. Release of the silencing effect of HNS on its promoter was previously demonstrated for Aar in EAEC 30 . We are currently characterizing the hierarchy of Rnr-interactions in the context of HNS and AraC/XylS regulators in space and time to ascertain the regulatory landscape of Rnr.
Interestingly, AraC/XylS and HNS are critical regulators of fimbriae and pili in most pathogens. These essential organelles allow bacteria to colonize their hosts and persist on abiotic surfaces. In EAEC, the AraC family member AggR controls the expression of AAF 23,24 . In ETEC, the expression of many types of adhesive pili www.nature.com/scientificreports/ depends on the AraC member Rns 39 . In contrast, ToxT, an AraC member from V. cholera, regulates the transcription of genes encoding the two major virulence factors, the toxin-coregulated pilus (TCP) and cholera toxin (CT) 40 . Unlike EAEC, EPEC relies on virulence factors that promote LA and intimate contact with the epithelium mediated by a T3SS, which results in attaching and effacing (AE) lesions, the hallmark of EPEC infection. The LA phenotype is determined by the BFP, which is regulated by PerA (a member of the AraC/XylS family) [41][42][43] . In this study, we show that Rnr controls the expression of PerA and the expression of the adhesion protein intimin and its receptor Tir (Fig. 4).
Most importantly, we found that both Aar and Rnr impact EAEC and EPEC colonization of human intestinal organoids (Figs. 8 and 9). Although ANR members may impinge on multiple regulons within the same bacterial cell to modulate multiple virulence traits, ANR's primary function could be perhaps the modulation of colonization dynamics of enteric pathogens in the intestine through repression and anti-repression of AraC/ XylS and HNS during pathogen-host interactions, and that this process may enable bacteria to spread out and persist in the host. Additional studies in suitable animal models should be conducted to test this hypothesis. In summary, our findings suggest that despite the low homology between ANR cognates in pathogens with different www.nature.com/scientificreports/ pathogenic schemes, ANR is a highly conserved mechanism of regulation of AraC/XylS and HNS regulators in diarrheagenic bacteria (Fig. 10).

Materials and methods
Bacterial strain and growth conditions. Bacterial strains and plasmids used in this study are listed in Table 1. 042aar was previously generated in our laboratory 29 . The rnr locus (2,721,409-2,721,642, GenBank FM180568.1) in EPEC E2348C_2635 was replaced with kanamycin (km) by lambda red recombination 44 . Our strategy, replaced the rnr locus with the Km gene from start to stop codon to avoid polar effects in downstram genes. Briefly, EPEC(pKD46) was grown in 100 ml of LB at 30 °C to an OD 600 ηm of 0.4. The lambda red enzyme was induced with 20 mM of l-arabinose for 1 h. The culture was centrifuged, and the bacterial pellet was prepared for electroporation using sterile water. 100 µl of electrocompetent cells were mixed with ~ 500 ng of DNA and electroporated using a Gene Pulser Xcell Bio-Rad electroporator. Cells were recovered in 3 ml of SOC medium at 37 °C for 3 h. The DNA used for the electroporation was amplified by PCR using EPEC as a template and specific rnr flanking primers (    45 . Counts for each annotated genomic feature were determined by htseqcount (http:// www-huber. embl. de/ users/ anders/ HTSeq/ doc/ count. html). Differential expression between counts for each feature was then calculated with DESeq 46 using the false-detection rate-adjusted Benjamini Hochberg P value. The P-value obtained from the test was corrected, and the false discovery rate (FDR) was used as a key indicator of differentially expressed genes. The fold change of differentially expressed genes vs. P value was plotted using GraphPad Prism 6 (GraphPad Software, Inc., CA, USA). During the analysis, fold change ≥ 2 and FDR < 0.05 were set as screening criteria. Fold change indicates the ratio of expression levels between the samples.

Real-time quantitative reverse transcription-PCR (qRT-PCR).
Overnight bacterial cultures of EPEC E2348/69 derivatives were diluted 1:50 in 13 ml of DMEM-HG (ANR-inducing conditions) and incu-  www.nature.com/scientificreports/ bated at 37 °C without shaking for 2, 3, and 4 h. Extraction of RNA, cDNA synthesis, and qRT-PCR assays were performed as previously described 23 . Primers for the qRT-PCR analysis are reported in Table 2. Reactions were run in experimental duplicate using two independent cDNA preparations. Expression levels for each queried gene were normalized to the constitutively expressed rrsB of EPEC as previously described 23,47 . Detection of AAF/II, EspA, and EspB. To detect the major fimbria subunit of AAF/II (AafA), strains were grown in 13 ml of DMEM-HG to reach an OD 600 ηm of 0.8. Bacteria were pelleted, resuspended in 100 µl of 0.5 mM Tris, 75 mM NaCl and heated for 30 min at 65 °C. AafA was analyzed in heat-prep supernatants by SDS-PAGE and Western blot analysis. Protein samples were separated in acrylamide gels and transferred to Immobilon-P membranes (BioRad, Hercules, CA, USA) using standard protocols. Membranes were incubated overnight with a polyclonal anti-AafA antibody. The next day, membranes were washed twice in PBS-0.1% tween and incubated for 1 h with horseradish peroxidase-conjugated goat anti-rabbit IgG antibody (ThermoFisher). Membranes were washed thrice with PBS-0.1% tween and imaged using ChemiDoc Imaging System (BioRad). Expression of EspA and EspB was analyzed in whole-cell and supernatants of EPEC cultures by Western blot. Briefly, EPEC derivatives were grown overnight in LB, diluted 1:50 in 13 ml of DMEM-HG, and grown at 37 °C to an OD 600 ηm optical density of 1.0. The cultures were centrifuged at 20,000×g for 5 min, and bacterial pellets were prepared for Western blot analysis. On the other hand, supernatants were collected and filtered through a 0.22 μm filter (Millipore). The supernatants were precipitated with 10% trichloroacetic acid overnight at 4 °C and centrifuged at 18,000×g for 15 min. Protein pellets were washed with 1 ml of acetone and dissolved in SDS-PAGE sample buffer saturated in Tris. Proteins were analyzed by Western blot using polyclonal anti-EspA and anti-EspB antibodies (kindly donated by Dr. James Kaper and Jane M Michalski, University of Maryland). Cytoplasmic GroEL (Abcam) was detected as an internal loading control and as an indicator of cytoplasmic protein contamination in the secreted protein fraction.
Human intestinal organoid culture. Human intestinal organoids used in this study were derived from a colonoid/enteroid repository previously established from deidentified biopsy specimens from healthy subjects who provided written informed consent at Johns Hopkins University by approved guidelines and regulations (IRB NA_00038329).
The maintenance of human organoids and preparation of colonoid monolayers were previously described 36 . Briefly, organoids were routinely cultured as 3D cysts embedded in Matrigel (Corning) and used to prepare cell monolayers in 24-well, 0.4 µm pore size polyester membrane cell culture inserts (Transwell supports, Corning) precoated with 100 µl of 34 µg/ml of human collagen IV solution (Sigma). Intestinal monolayers were routinely grown at 37 °C with 5% CO 2 until confluency as assessed by the increase in transepithelial electrical resistance (TEER), measured using an epithelial volt/ohm meter (EVOM, World Precision Instruments). Confluent monolayers were differentiated for five days before infections.
Confocal microscopy. EAEC 042 and EPEC E2348/69 derivatives were grown overnight in LB supplemented with appropriate antibiotics (Sigma Chemical Co, St. Louis, MO). The next day, overnight cultures were diluted 1:50 (V/V) in DMEM-HG medium (Invitrogen, USA) and incubated at 37 °C with shaking to the midlog phase (OD 600 = 0.6) to induce ANR expression. Bacterial cultures were adjusted to 10 8 CFU/ml in PBS, and 10 µl (10 6 CFU) was added to the apical surface of colonoid monolayers. Cells were infected for six hours at 37 °C, 5% of CO 2 . Following the bacterial infection, the cells were fixed with Carnoy's solution (90% methanol, 10% glacial acetic acid), washed three times with PBS, permeabilized with 0.1% saponin, and blocked with 2% bovine serum albumin/fetal bovine serum for 30 min (Sigma Aldrich, USA). Cells were rinsed with PBS and incubated overnight at 4 °C with primary antibodies diluted 1:100 in PBS containing 15% FBS and 2% BSA. As primary antibodies, we used anti-O44 for EAEC staining, anti-intimin for EPEC staining, and Alexa-phalloidin for cell-actin staining. Stained cells were washed three times with PBS, followed by incubation with appropriate Alexa-conjugated secondary antibodies (Molecular Probes/Invitrogen, USA) diluted 1:500 in PBS. Hoechst (Vector Laboratories, USA) was used at a 1:1000 dilution in PBS for nucleus/DNA labeling. After incubation, cells were washed three times for 5 min each and mounted in ProLong Gold (Vector Laboratories, USA) over- www.nature.com/scientificreports/ night at 4 °C. Confocal imaging was carried out at the Imaging Core Facility at the University of Virginia using an LSM-710 Multiphoton laser-scanning confocal microscope (Zeiss, Germany) running ZEN 2012 (black and blue edition) imaging software (Zeiss, Germany). Images were captured with a 64 × oil objective. The same settings were used to image across samples for quantitative analysis of bacteria and AE lesions (actin polymerization). At least 10 Z-stack images of 2-μm intervals up to a depth of 20 μm for each treatment were obtained from three independent experiments. The Z-stack image depth corresponded to a detectable green fluorescent signal (bacteria) on colonoids stained for actin. A robust fluorescent signal (polymerized actin) corresponds to AE lesion formation in EPEC-infected bacteria. The relative number of bacteria (or AE lesions) was measured in randomly selected Z-stack microscopic fields by ImageJ software (NIH), using the particle enumeration algorithm as previously reported 36 . Briefly, single images from Z-stacks exhibiting green fluorescent bacteria or fluorescent polymerized actin (AE lesions) were obtained using Blue edition ZEN2012 software Zeiss, (Zeiss, Germany). Images were opened as 16-bit type images with ImageJ. Threshold values were adjusted to eliminate the background. Bacteria or AE lesions (particles) were enumerated in images processed as Binary > Watershed images. This algorithm separates particles that are close together (e.g., aggregated bacteria). Lastly, images were analyzed as particles set as size (pixel2) = I0-infinite, which is relatively close to the size of E. coli in 64 × confocal images. Particle (bacteria or AE lesions) counts in each image of a Z-stack were input in an Excel sheet and plotted using the Prism software (Graph Pad).
Bioinformatic and statistical analysis. The amino acid homology of ANR proteins was obtained using Clustalw algorithms (http:// www. genome. jp/ tools/ clust alw/). Values are presented as mean ± and standard error of the mean (SEM). Statistical significance was determined using analysis of variance (ANOVA) with Bonferroni's post-test (Prism GraphPad) to compare groups, including a minimum of n = 3 replicates. Results were considered significant at P < 0.05.

Institutional guidelines.
Human intestinal organoids used in this study were derived from a colonoid/ enteroid repository previously established from deidentified biopsy specimens from healthy subjects, who provided written informed consent at Johns Hopkins University (Protocol NA_00038329). All methods were carried out in accordance with the University of Virginia-approved guidelines and regulations (IRB-HSR # 18959). All experimental protocols were approved by the University of Virginia Review Board (IBC number: 1329-11).

Data availability
The RNA-seq data discussed in this publication have been deposited in NCBI's Gene Expression Omnibus and are accessible through GEO accession number: GSE225626. To review GEO accession GSE225626: go to https:// www. ncbi. nlm. nih. gov/ geo/ query/ acc. cgi? acc= GSE22 5626. Enter token yhczigeijvsnnir into the box.